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Abstract 

The HERA experiments HI and ZEUS have recently presented a number of interesting results 
from studies of the hadronic final state produced in ep collisions. These studies have focussed on 
deep inelastic scattering events and on the photoproduction of jets, and also on the comparison of 
the final states produced in these two regimes. Some of these results are reviewed here. 
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1 Introduction 


Analyses of the distribution of final state hadron^ in ep scattering events can provide tests of the theory 
of strong interactions, quantum chromo-dymamics (QCD). These studies have yielded information on 
the value and running of ag, on the physics of parton shower evolution, on the structure of the photon 
and proton and on the underlying parton scattering dynamics. There are too many recent results 
from HI and ZEUS to attempt a complete review here. Instead I will concentrate on the measurement 
of fragmentation functions in deep inelastic scattering (DIS) Q, some results from the measurement 
of jet cross sections in photoproduction ||^ |3|, a study of prompt photon production and some 
results from the comparison of the hadronic final state produced in deep inelastic scattering with 
that produced in photoproduction i i- These results pertain to data collected in the 1994 and 
1995 periods of data taking, during which HERA provided collisions of 27.5 GeV positrons against 
820.0 GeV protons, delivering about 3 nb“^ and 6 nb~^ of luminosity respectively. 

2 Hadronic Final State in Deep Inelastic Scattering 

In the leading order quark parton model of DIS, the incident positron scatters from a quark in the 
proton via the exchange of a boson whose negative invariant mass squared is denoted Q^. The quark 
carries a fraction XBj of the proton’s longitudinal momentum and energy into the hard subprocess. 
This simple picture is sufficient to allow the appreciation of the results presented here. 

2.1 Fragmentation Functions 

To characterize the parton showering and subsequent hadronization which determine the distribution 
of hadrons in the final state, the distribution of the scaled momentum Xp is measured. The scaled 
momentum is defined as the magnitude of a hadron’s momentum relative to the maximum momentum 
available to it, Xp = Phadron/pmax- The measurement of scaled momentum distributions is carried 
out in the “Breit” or “brickwall” frame where pmax = — <5/2. In this frame the current hemisphere in 
DIS is directly comparable to one hemisphere in an e'’'e“ annihilation event where the centre-of-mass 
energy E* would represent the energy scale, Q, and pmax = E*I2. 

In fact a variable transformation is commonly made to ^ = ln(l/xp) which has an approximately 
Gaussian distribution. The distribution of ^ is called a fragmentation function. As an example, 
fragmentation functions in two bins of are shown in Figure for data taken by HI during the 
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Figure 1: Fragmentation functions for (a) 12 < < 80 GeV^ and (b) > 100 GeV^. Gaussian fits 

to the distributions are superimposed. 

1993 running period 0]. The data support the intuitive expectation that as increases, the average 
charged multiplicity should increase, {xp) should decrease, and (^) should increase. 

To quantify the evolution of ln(l/xp), a Gaussian fit of the ^ distribution has been made for 
several different Q bins. The resulting fit parameters for the peaks of the Gaussian distributions, ^peakj 


^In the ep world it is common to refer to ail of the final state particles other than the scattered lepton as the “hadrons”. 




are shown in Figure [||. The data are consistent with a QCD inspired picture of fragmentation in 
which ^peak IS expected to scale approximately linearly with InQ. A pure phase space model would 
entail a slope of unity in the dependence of ^peak on InQ while a fit to the data yield 0.75 ± 0.05. 
Therefore the data prefer coherent fragmentation as predicted by QCD. (See also reference and 
references therein.) 
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Figure 2: Peak values of the Gaussian fit to the ^ = ln(l/xp) distribution versus Q for HI 1993 (open 
squares) and HI 1994 (filled squares) data. The line shows a linear fit to the data. Also shown is the 
corresponding ^ distribution parameter versus E* for e'’'e“ data (circles). 

Results for the evolution of the ^peak fit parameter with E* for e"*“e“ data are also shown in 
Figure The consistency of the e'’'e“ data with the ep data supports the application of a single 
universal treatment of parton showering and hadronization. 

3 Photoproduction of Jets 

Deep inelastic scattering processes (Figure ^(a)) are characterized by rather high values of the negative 
squared four momentum of the boson (usually a photon), > 0(1) GeV^. In photoproduction on 





Figure 3: Diagrams for DIS (a) and direct (b) and resolved (c) photoproduction. 

the other hand the photon is nearly real. When jets are produced in photoproduction, their transverse 
energy sets the energy scale of the process, also often represented by Q^. 

At leading order there are two photoproduction processes which lead to the production of high 
transverse energy jets. In “direct” photoproduction the photon participates directly in a hard scat¬ 
tering with a parton from the proton as illustrated in Figure ^(b). In “resolved” photoproduction 
(Figure ^(c)) the photon first fluctuates into a hadronic state a parton from which enters the hard 
subprocess. It is therefore necessary to introduce another Bjorken scaling variable, denoted x^, which 
represents the fraction of the photon’s momentum which is carried by the interacting parton. 







3.1 Jet Cross sections 


ZEUS has measured the dijet photoproduction cross section differential in the average jet pseudora¬ 
pidity, da/dfj, for events with at least two > 6 GeV jets [^. By making the requirement that 

> 0.75, we effectively examine the direct photoproduction cross section |^]. The results are shown 
in Figure ^(a) where three different jet finders have been used to define the cross section. EUCELL 
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Figure 4: Dijet photoproduction cross section, da/dfj, for (a) > 0.75 and (b) 0.30 < x^ < 0.75. The 

different symbols show the measured cross section using three different jet finding algorithms. The 
curves show NLO calculations of the cross section using the CTEQ3M proton parton distributions 
and the GRV photon parton distributions for two different values of the Rsep parameter (see text). 
(The shaded band shows the correlated uncertainty on the absolute energy scale of the detector.) 


and PUCELL are cone algorithms which satisfy the Snowmass convention |Q. The difference between 
them illustrates the sensitivity of the measurements to the choice of seed particle and the treatment 
of overlapping cones. KTCLUS is a cluster jet finding algorithm |T^ . 

The results are compared to next to leading order (NLO) predictions for the cross sections using 
the KTCLUS algorithm and two different values of the Rsep parameter [11|. (Rsep defines the farthest 
distance in pseudorapidity and azimuth that two partons can have from one another and still be 
combined into one jet ||l^.) The prediction shown by the dashed line may be directly compared with 
either the PUCELL or KTCLUS cross sections which are equivalent at NLO and have an effective 
Rsep value of 1. The theory provides a reasonably good description of both the normalization and 
shape of the measured distribution. In EUCELL the effective value of Rsep is not precisely defined 
but is in the range 1.5—s■2. The theory is therefore able to describe the spread of the cross sections 
due to the different jet finding algorithms. 

The situation is considerably worse for the resolved photon processes. The cross section da/dfj 
for 0.30 < x^ < 0.75 is shown in Figure |^(b). Here the NLO calculation for the KTCLUS algorithm 
(dashed line) does not describe the measured KTCLUS cross section. Also the cross section suffers 
more variation under change of jet algorithm and this variation is not describable by the NLO calcula¬ 
tion. However it must be noted that much of the uncertainty due to the photon structure is contained 
in the resolved sample and that also the higher order processes which are not included in the fixed 
order calculation would show up at low x^. 

It is possible to define an effective parton distribution of the photon, = X^(g + g) + 9/4g, which 
can be measured in photoproduction |13|. Note that in contrast to deep inelastic ey experiments which 
only determine the quark content of the photon, the HERA jp measurements are sensitive also to the 
gluonic term of fi/y. By comparison of measured parton level cross sections with the PYTHIA |14] 
Monte Carlo leading order prediction HI have extracted the effective parton distribution of the photon 






















in three bins of as a function of the hard scattering scale, (the square of the partonic transverse 
momentum). |^] This is shown in Figure 
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Figure 5: The symbols show the photon parton density versus in three bins of x^. The curves 
represent the photon parton density calculated using the GRV-LO parton distribution functions. 


The data are consistent with a logarithmic rise of with We recall that a calculation of 

the quark content of the photon from the box diagram 77 * ^ qq yields the form = a{xry) + 

b{x^)logpf 1^]. It is therefore interesting to speculate whether the logarithmic rise of fi/^ with p'^ 
reveals a component of the resolved photon which results directly from a splitting of the photon into 
a qq pair. This has been called the “anomalous” component in contrast to that part of the photon’s 
structure which behaves as though the photon has fluctuated into a vector meson. 


3.2 Prompt Photon Production 

A particularly interesting photoproduction event topology occurs when a high transverse energy pho¬ 
ton takes the place of one of the jets. Such processes are called prompt photon processes. These events 
provide a clean QCD laboratory for the study of quark jet showering and fragmentation and for the 
measurement of the quark content of the resolved photon. This is also the best process for the study 
of initial-final state QCD coherence. In direct photoproduction there is a single flow of colour from 
the incoming to the outgoing quark, in contrast to the situation in generic dijet production where 
many different possible colour flows contribute to each event topology. 

To select a sample of prompt photon events a jet and an isolated electromagnetic energy deposit 
in the ZEUS barrel calorimeter with no associated charged particle found in the tracking chamber, 
are required. However some dijet events, where one jet consists largely of a single vr^ or rj, can pass 
this selection. This contribution can be understood using a Monte Carlo simulation as shown in 
Figure ^(a) Q. The fraction of the electromagnetic shower energy which is contained in the highest 
energy calorimeter cell is denoted /max- Photons give signals with higher values of /max than the vr® 
background and it is possible to describe the measured distribution of /max by allowing for a small 
component. This also allows for an admixture of r/’s. 

The x^ distribution for prompt photon events after subtraction of the vr^ background is shown in 
Figure ^b). The distribution is roughly consistent in shape and normalization with a leading order 
QCD prediction. Most of the events have high values of x^ and are therefore associated with direct 
prompt photon production. We look forward to a larger sample of prompt photon events with the 
1996 and 1997 delivered luminosities. 
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Figure 6: (a) The distribution of /max for the data (solid histogram), the simulated prompt photon 
events (dashed histogram) and the simulated background (dotted histogram). In (b) the black 
dots show the distribution for prompt photon events after background subtraction. The dashed 
histogram shows the expectation from the PYTHIA Monte Carlo simulation. The dotted histogram 
shows the simulated distribution separately for the (leading order) resolved component. 


4 Comparison of Deep Inelastic Scattering and Photoproduction 

Two interesting recent results compare the hadronic final states in deep inelastic scattering and pho¬ 
toproduction. These studies are particularly appealing to the experimentalist since it is possible to 
make some interpretation of the data which does not rely on a comparison of the data to theory. 

4.1 Transverse Energy Flow 

HI has measured the flow of transverse energy in the hadronic centre of mass frame across a large 
range of Q^, as shown in Figure ^ Q. The hadronic centre of mass frame pseudorapidity, r]*, is defined 
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Figure 7: Average transverse energy, Et, in bins of the hadronic centre of mass frame pseudorapidity, 
1 ]*, and (Q^). 

such that the proton direction is negative so the current region is at positive r]*. 

The Et in the current region of course grows with (Q^). However the remarkable feature of 
this distribution is that the mean transverse energy in the central region, r/* ~ 0, is practically 
independent of from the soft photoproduction region, (Q^) = 0, up to the very deep inelastic 
region, (Q^) = 520 GeV^. HI therefore find that there is a smooth transition in the development of 
the global flow of transverse energy between the photoproduction and DIS regimes. 




























4.2 Jet Shapes 

A complementary result from ZEUS concerns the internal structure of jets Q. The jets are found using 
a cone algorithm with jet radius R = 1. The jet shape parameter, has been defined as the fraction 
of the jets transverse energy which is contained within a smaller concentric cone of radius r. This is 
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Figure 8: Jet shapes in dijet photoproduction for high (open circles) and low (black dots) samples. 
The PYTHIA predictions for the low and high regions are shown by the solid and dashed lines 
respectively. Thicker lines show the PYTHIA predictions where multiple interactions between the 
constituents of the photon and proton have been allowed. 

shown in Figure ^ for a sample of photoproduction events which have two jets of > 14 GeV 

and — 1 < < 2. The jet shape is shown separately for direct and resolved photon events. The 

jet shape for the x<y > 0.75 sample is narrower than the shape of the jets for the resolved photon 
sample. This effect is reproduced by the parton shower Monte Carlo model PYTHIA. Allowing for 
secondary interactions between the photon and proton remnant particles in the simulation improves 
the description of the x^ < 0.75 data. 

One can speculate that the observed difference in the high and low x^ jet shapes is due primarily 
to the colour charge of the parton which initiates the jet. In direct photoproduction this parton is 
predominantly a quark whilst in resolved photoproduction the jets are much more often initiated by 
gluons. If the direct photoproduction jets are narrower than the resolved photoproduction jets because 
of the predominance of quark jets, then one would expect the jet shapes of the high x-y events to be 
comparable to jet shapes in DIS. This is indeed the case as shown in Figure |^]. The high x^ jet 
shape is shown in four bins of and is found to be comparable with the jet shape for DIS events 
(Q^ > 4 GeV^) having at least one Et > 14 GeV jet p. 


5 Conclusions 

Quantum chromodynamics is a successful theory. For a number of years it has been clearly established 
that in the perturbative regime QCD matrix elements describe the dynamics of parton scattering. 
Here it is shown that ep scattering provides the means to perform detailed tests of the QCD inspired 
models for parton showering and hadronization. These models, as implemented in Monte Carlo event 
generators, are able to describe the measured observables. There is also strong evidence that the 
fragmentation process has a universal nature given the similarity of the hadronic final states in e+e“ 
and in DIS events as well as in DIS and photoproduction events. 
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Figure 9: Jet shape in direct photoproduction (dark dots) compared with jet shape in DIS (light 
dots) in four bins of . The data are not corrected for detector effects and are compared with 
PYTHIA Monte Carlo predictions to which a full simulation of the ZEUS detector has been applied. 
The prediction for the the high jet shape is shown by the dark line and the prediction for the pure 
leading order direct processes is shown by the light line. 
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